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Abstract. - We show that the periodic modulation of the Hamiltonian parameters for ID cor- 
related fermionic systems can be used to parametrically amplify their bosonic collective modes. 
Treating the problem within the Luttinger liquid picture, we show how charge and spin density 
waves with different momenta are simultaneously amplified. We discuss the implementation of our 
predictions for cold atoms in ID modulated optical lattices, showing that the fermionic momentum 
distribution directly provides a clear signature of spin-charge separation. 



^ Introduction. — The swing is the best known exam- 
^ple of a classical system showing parametric resonances 
. The periodic modulation of the effective swing length 
induced by the motion of legs leads to the exponential am- 
^ plification of the oscillations if the modulation frequency is 
^ chosen commensurately with the natural frequency of the 
swing. Quantizing this classical problem as a harmonic os- 
(— I cillator with modulated parabolic confinement leads to the 
O appearance of an exponential divergence in the time evo- 
^ O l ution of the bosonic rising and lowering operators. This 
effect is particularly strong if the modulation is around 
P\| twice the natural oscillator frequency. 
> In a modulated system of many bosonic oscillators only 
those fulfilling the resonance condition will be ampli- 
. fied, making parametric resonance a spectroscopic tool in 
many-body quantum systems. These ideas acquired par- 
^^ticular relevance since cold atoms in optical lattices have 
been realized [2] . As the intensity of the lattice can be fully 
controlled by the laser power one can study parametric 
^^modulations in correlated quantum systems with current 
^experimental tools [s]. Similarly, the periodic modula- 
" tion of the transverse confinement in cigar-shaped Bose- 
r> Einstein condensates has been shown to induce the para- 
^ metric amplification of Faraday waves [i] . From the the- 
oretical point of view, parametric amplification of Bogoli- 
ubov quasiparticles for bosonic clouds in optical lattices 
have been already investigated in the past (5|[6]. In these 
systems, the bosonic nature of quasiparticles appears al- 
ready when interactions are treated at mean-field level, 
and allows for the amplification to occur. 

In this paper we analyze parametric resonances in many- 



body fermionic systems, starting with the very question 
whether the amplification can occur at all. Indeed, in 
contrast to the bosonic case, in fermionic systems any 
mean-field treatment of interactions, including the pres- 
ence of broken symmetries, preserves the fermionic nature 
of quasiparticles. The Fault principle thus blocks their am- 
plification, as can be easily checked by direct calculation 
. Then the question rises if bosonic collective excitations 
of a fermionic many-body system can be subject to am- 
plification by modulating a parameter in the microscopic 
Hamiltonian. In order to address this fundamental ques- 
tion we need to treat correlations in a fermionic system 
beyond mean-field level. In this work we thus confine our 
investigation to one-dimensional (ID) correlated fermions 
within the Luttinger liquid picture, in which interactions 
are treated exactly and the system is naturally diagonal- 
ized in terms of collective bosonic spin and charge density 
waves [7][8] . According to the Luttinger liquid theory these 
modes disperse with two different group velocities, giving 
rise to the so-called spin-charge separation (see fig. [T]). 
This fundamental issue in condensed matter physics has 
been detected in transport experiments on quantum wires 
[oj and by angle-resolved photoemission spectroscopy of 
ID SrCu02 |10 only recently. Due to their great tunabil- 
ity, cold atomic gases in optical lattices are also promising 
candidates for the experimental detection of spin-charge 
separation in ID systems, as shown by several theoretical 
proposals 
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Here we show how a spatially homogeneous time- 
periodic modulation of the intensity of the optical lattice 
indeed leads to the amplification of charge and spin den- 
sity waves of a ID correlated cloud of ultracold fermionic 
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Fig. 1: Sketch of the dispersions ujq of the charge and spin den- 
sity waves as a function of momentum q, with group velocities 
v'^ and v'^ , respectively. The pumping frequency Q. amplifies col- 
lective spin and charge modes in the vicinity of the momenta 
qc and q^. 



atoms. If 51 is the modulation frequency, the charge and 
spin waves with energy v'^q^ — v'^q^ = 51/2 will be ampli- 
fied (see fig.[l]). Due to the different group velocities v'^ and 
v'^ for the charge and spin channels, respectively, the reso- 
nant condition amplifies different wavenumbers for the two 
branches, qc and qs- On top of showing the feasibility of 
parametric amplification in correlated fermionic systems, 
we also propose this technique as a tool to systematically 
investigate the spin-charge separation in experiments. In- 
deed, we discuss the effect of the amplification above on 
the fermionic momentum distribution and show how the 
latter exhibits well defined shoulders directly related to 
the wavenumbers Qc and qs- These structures are partic- 
ularly evident after not too long modulation times. As 
the momentum distribution is the standard quantity mea- 
sured in time-of-flight experiments on cold atomic clouds 
[2], our analysis has direct implications on investigations 
of correlated fermionic systems with current experimental 
tools. 

Model. — We consider interacting fermionic atoms 
with (pseudo-)spin 1/2 confined into ID cigars (as, e.g., 
realized in optical lattices [2]), with a further periodic po- 
tential along the ID axis. Despite this additional poten- 
tial, we assume the atoms to be in their metallic phase, as 
opposed to the recently investigated Mott-insulator phase 
18 -20 . In order to obtain analytical results, we disregard 



trapping and finite size effects, as they will not qualita- 
tively modify our results. The parametric excitation of 
the system is achieved by periodically modulating the in- 
tensity of the optical lattice along the ID system, thereby 
shrinking the Wannier wavefunctions associated to each 
lattice site. This has the two-fold effect of modulating the 
hopping rate between neighboring sites (i.e., the kinetic 
energy) as well as the on-site repulsion for multiple oc- 
cupancy. As the former is exponentially sensitive to the 
wavefunction overlap it is the dominant parametric mod- 
ulation term in the problem and we will focus on this for 
simplicity. It should be however noticed that the inclusion 
of the smaller parametric modulation of the interactions 
would not significantly affect the consequences discussed 
in the following. 



Once translated into the Bloch-band language, and if 
we focus on the low-energy physics of the problem, the 
modulation in the kinetic term yields an effectively time- 
dependent Fermi velocity for the atoms close to the Fermi 
level, VF{t) = v-p + SvF{t). Here Sv-p(t) = 7tiF sin (51i) , 
where v-p is the Fermi velocity at equilibrium (i.e., for 
times t < before the modulation starts) , and 51 and 7 are 
the frequency and intensity of the parametric modulation, 
respectively. Our focus on the low-energy sector of the 
many-body problem naturally suggests a Luttinger liquid 
approach [Tjjs] to the correlated system, which allows for 
an exact treatment of interactions. We thus linearize the 
single-particle spectrum in the vicinity of the Fermi level 
for momenta ||fc| — fepj < A, with kp the Fermi wavevector 
and A an ultraviolet cutoff. The time-dependent Hamil- 
tonian describing the system of size L therefore readf[^ 
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with fkr = rk — kp. Here, c^^^ {ckar) creates (an- 
nihilates) a T-moving fermion with momentum k and 
spin (7 =t,i (r = 1 (—1) corresponds to right (left) 
movers) while pj"' = "^k^k-q a T'^k,(j.T is the correspond- 
ing fermionic density operator. In ([T]) we keep the general 
form of the interaction between fermionic densities with 
generic spin and branch indices. Within the standard Lut- 



tinger liquid theory we have V., 
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V?%,q = - Vi^, V:;i:^^ = Fo^, where Vlj'^ is the 
Fourier transform of the microscopic interaction between 
fermions with parallel/antiparallel spins. This allows to 
treat short range interactions (like contact s-wave scat- 
tering for neutral fermions) where Pauli principle imposes 
Vq = 0, as well as finite range spin-invariant ones {e.g., 
between charged fermions) with Vq = V^. Eq. ([T]) does 
not include backscattering between particles with oppo- 
site spins as well as umklapp scattering as those terms 
are usually negligible at equilibrium and away from half- 
filling [7][8] . The effect of the smaller periodic modulation 
of these interaction terms on the energy absorption has 
been considered in ref. [2l] . 

Introducing the bosonic operators for charge and spin 
density fluctuations 6^ = (7r/L|g|)i/2 e(rg)(p^t + pj^) 

and h\ = (7r/-^-|g|)^/^ Er " Pg^)' respectively, 

eq. ([1]) transforms into the separable Hamiltonian [t] 

H{t) = ^ \q\ [Al{t)bfhl + B- {hfh% + 6^,5^)" , 

a— c,s 

(2) 

with Al{t) = vp{t) + (yj' + yo^)/2^, Al{t) = vp{t) + 
(T/ll - Fo^)/27r, = (l/ll + V,^ - K,i,)/4^- and B^ = 
(yii - - 4j/4^. 

^In the remaining of the paper, we set h = k-Q = 1. 
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The time- independent part in ([2| can be diagonahzed by 
means of the Bogohubov transformation &° — cosh /3° -|- 
sinhc/?^ in terms of the bosonic fields /3°, such that 
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c.,«+<5«F(i)|g|cosh(2¥.«)) /3f A" 



q\ sm 



h(2^^)(/3f/3^1, + /3!!/; 
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with 



velocities are v" 



The different charge and spin group 
g - [A^{Of - ABf]^/'^. For short range 
interactions, one can neglect the momentum dependence 
of the group velocities such that ~ v"", and the dis- 
persion of the spin and charge density waves is linear 
(see fig. [l} . The coefRcients of the Bogoliubov transfor- 
mation read sinh(/75 = -[yl^(0)/2t;^ - 1/2] ^/^^ coshv?;^ = 
[A;;(0)/2w^ + l/2]i/2^ sinh(^^ = [Al{0)/2v^g - 1/2]^/^, and 
coshv?^ = [Al{0)/2vl + l/2]i/2. 

The crucial point to notice at this level is that the para- 
metric modulation introduces a time-dependent anoma- 
lous term in the Hamiltonian ([S]), creating and annihilat- 
ing pairs of bosonic charge and spin density waves. As 
the modulation is homogeneous in space (i.e., at zero 
wavenumber), the new terms create or annihilate pairs 
of excitations with opposite wavenumber, as requested 
by momentum conservation [6j. In addition, the induced 
anomalous terms are proportional to sinh(2iy9^) and thus 

correctly vanish in the non- interacting limit vj'''^ = 0, 
where fermionicity forbids parametric amplification. 

Parametric amplification. From the Hamilto- 
nian ([3| above we can now determine the time evolution 
of the operators /3°, showing that indeed parametric 
amplification of collective bosonic modes in a fermionic 
system is possible. With eq. ([s]), the Heisenberg equation 
of motion for the operator reads $g{t) — — i[Wq + 
SvFit)\q\coshi2^-gWg{t) - i(5«F(t)|<z|sinh(2(p«)/3^;(i). 

Defining /3^{t) = f3^(t), 
assuming a weak parametric modulation (7^1) and 
retaining only slow terms near the resonance (rotating 
wave approximation, i.e., for in the vicinity of 2a;°), 
the_equation of motion and its adjoint transform into 

0, where 




■,q — 7Vf|<z| sinh(2i^^)/2. Solving this equation with 
the appropriate initial conditions (3^(0) — /3q(0) and 
/3^(0) - 



^ql3l\{0), we obtain 
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where a;^ = - fi/2 ± 



we defined 



;--n/2y-^f. In Q, 
and = ^^/{uj-g_ 



i^q^)- Thus, in a narrow "resonant window" of energy 
\ujg — il/2| < l^^l, the frequencies ujq± acquire an imagi- 
nary part, leading to the exponential amplification of the 



corresponding bosonic modes. Outside this window, the 
modes evolve according to their coherent dynamics and 
are therefore not amplified. Indeed, out of eq. Q, it is 
easy to verify that, on- and off-resonance, the evolution of 
the Bogoliubov operators is 



f3qit) 
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sinh(|C|t)/3^t,(0) 
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Fermionic momentum distribution. The results 
above prove the possibility of amplifying bosonic collective 
modes (with different wavenumbers for charge and spin 
modes) while parametrically modulating the underlying 
ID fermionic many-body Hamiltonian. However, the de- 
tection of this amplification is not necessarily easy from 
the experimental point of view. In the case of cold atoms 
in optical lattices this would require a (spin-resolved) mea- 
surement of the cloud density during the parametric mod- 
ulation, without opening the trap. Very recently, in-situ 
measurements on confined ultracold atomic gases have 
been reported [ij 22 23] . The spatial resolution of these 
measurements would allow for the detection of density 
modulations in the cloud associated to the parametric am- 
plification of charge density waves. The detection of spin- 
charge separation, in addition to the rich physics of the 
Hubbard model for cold atoms in optical lattices, could 
motivate further experimental efforts towards the realiza- 
tion of local spin-resolved in-situ measurements not re- 
ported so far. 

In view of this experimental challenge we now discuss 
the consequences of our analysis on the fermionic momen- 
tum distribution, which is the standard quantity measured 
in time-of-flight experiments [2]. Due to the different am- 
plified momenta for charge and spin modes, we show that 
the fermionic momentum distribution shows clear signa- 
tures of spin-charge separation. 

The momentum distribution function (per spin 
channel) for r- movers is defined as nkrit) = 
f dxe^'^^ (ipl„{x,t)^rcr{0,t)) where, within our per- 
turbative treatment in 7 <C 1, (■ ■ • ) represents a thermal 
average with respect to the time-independent part of 
the Hamiltonian ([T]), i.e., for 7 = 0. The fermionic field 
operators creating a r-mover with spin a at position 
X and time t can be expressed in terms of the bosonic 
operators as [t] 



Vi(l - 



Ura{t), 



(6) 



with Ora{x,t) ^ [0';ix,t) + a9l{x,t)]/V2 (cr = I and 
— 1 correspond to spin up and down, respectively) and 
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the unitary Klein operator Ura increases the number 
of T-movers with spin a by one. 

With ([6|, we find after a lengthy but straightforward 
calculation 



L(l-C-2ViA)' 

Here, within the rotating wave approximation, we have 



(7) 



a,q>0 



+ [cos{qx) - 1] cosh (2^^) {{f3;it),pf{t)}) }, 

(8) 

with (3q{t) given in eq. Q. The summation over g in ([s]) 
is simplified by taking the evolution of the Bogoliubov op- 
erators in eq. ([s]) for wavenumbers belonging to the off- 
resonance and the narrow on-resonance windows. The 
latter are centered around the two resonant wavenumbers 
Qc/s = ^/2v'^/^ for the charge and spin modes, respectively. 

Our procedure allows for a fully analytical treatment of 
the problem, leading to the fermionic correlator 

i)^r.(0, t)) = {'lPl„{x),PraiO))o Aix, t) , (9) 

where the correlator without parametric amplification 
{i.e., for 7 = 0) in the zero temperature limit is [t] 



le 



27r(ra; + i0+) \x^ + 



A2 



(10) 



with a = (sinh ip'^ + sinh if^)/2. To obtain (10 1, we ap- 
proximated by its q limit ip"". This is justified pro- 
vided the sum over momenta g in ([s]) is cutoff at g ~ 1/A, 
where A is the screening length associated with the specific 
form of the interaction between the particles. This leads 
to the fermionic momentum distribution in the absence 
of the parametric amplification n1 ^. As the problem is 
fully symmetric between right and left movers, we focus 
on the former without loss of generality. For k > kp, the 
momentum distribution can be expressed as 



21/2-Q 



dQ'(Q + Q')""'^'K„_i/2(Q + Q'), 

(11) 

with Q = A(fc — fcp), r(z) and Ky(z) being the gamma and 
the modified Bessel functions, respectively. In particular, 
"fc-fep^ = 1 - "-fe+feF,R- The function is presented 
in fig. E] at time t = for contact and finite-range interac- 
tions. The latter may be of relevance for the treatment of 
trapped cold ions. 

The amplification factor in eq. ^ reads 



A{x,t) 



exp I E /i°(t)[cos(9„x)-l] j, (12) 
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Fig. 2; Fermionic momentum distribution for right movers 



nk,R{t) in the short-time regime (see eq. 15 I as a function of 
momentum k, scaled by the screening length A. In the fig- 
ure, T = 0, 7 = 0.1 and Xq^ = 1. (a) Case of contact s-wave 
scattering interaction, with Vq' /v-p = 0.5 and Vq — Vj^.^ — 0, 
such that Age = XqaVs/vc — 0.85 (b) Case of finite-range in- 
teractions, with Vq/vf = Vo^/up = 10 and Vji-p/wp = 2, such 
that A^c = XqsVs/vc — 0.31. Thicker fines correspond to larger 
times t. 

where 

h^it) = [1 + 2nB(r2/2)]K" cosh {2ip'') [cosh (k^^) - 1] , 

(13) 

and = 7|sinh(2v5'')|uF/2w''. Here, nB(a;) = (e"/^ - 
1)~^ is the Bose distribution at temperature T. It is im- 
portant to notice that the amplification factor equals 1 
in the non-interacting limit (where sinh((p'*) = 0). Once 
again, this highlights the importance of fermionic interac- 
tions as a necessary ingredient for the parametric amplifi- 
cation to occur. Out of eq. ^ the momentum distribution 
at finite times thus results in the convolution 
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dq 
2tt 



(14) 



with A{q,t) the Fourier transform of (12 1. 

Two qualitatively different regimes occur in the small- 
or large-time regimes, i.e., if h°-{t) <C 1 or h°'{t) ^ 1 in 
cq. ( |l2| ). In the first case, valid up to times of order ~ 
- In ([1 + 2riB(f^/2)]K° cosh {2ip'')) / K°-n, the expansion of 
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(12]) yields A{q,t) = 2TTS{q)- 
Qa) ^ 2(5((7)], leading to 



nk,K{t) 



'fc.R 



(15) 

Despite parametric amplification, the fermionic momen- 
tum distribution fulfills nk-kY,K = 1 — '^--fe+feF.R ^ i"^ the 
unperturbed case, guaranteeing particle-number conser- 
vation. For A: > /cp it shows two steps of size h'^/^{t)/2 
involving fermionic states with momenta up to q^/s away 
from the Fermi level, as exemplified in fig. [2] 

These are direct signatures of the parametric amplifica- 
tion of the charge and spin density waves and their obser- 
vation can thus be used to detect spin-charge separation 
in interacting ID Fermi systems. By spanning the exter- 
nal modulation frequency fi, the whole dispersion of the 
collective modes can be mapped. From the point of view 
of the measurement, our result is best visible in the "short- 
time regime" where the expansion above holds, leading to 
two well resolved steps of size up to order 1/2. This fact is 
crucial in order to experimentally detect the amplification 
in the momentum distribution against other smoothening 
factors, like, e.g.^ trapping and finite temperatures. 

Our analytical treatment allows formally the analysis 
of the "large-time regime" as well, where h°'{t) ^1. In 
this case the amplification factor ([12]) can be approx- 
imated as A{x,t) — A^{x,t)A^{x,t), with A°-{x,t) = 
J2n=-oci ^^P ( ~ (laX — 27rn)^/i°(t)/2). As a consequence, 
the Fourier transform results in 
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with f^it) = (l/y^lHt)) exp ( - mV2/i°(t)) leading to 



nk,K{t) = 2^ z rii 
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Thus, in the large-time limit the fermionic momentum dis- 
tribution shows many small-size steps stemming from both 
the charge and the spin sectors (see fig. [3|. Indeed, the 
form of /m(i) shows how for large h°-{t) more and more 
peaks of A{qjt) in (16) become relevant. This will limit 
the experimental resolution of the structures in nfc_R(t) 
in contrast to the short timescales. Moreover, at large 
times the exponential amplification of bosonic modes re- 
quires a treatment of their residual interaction beyond 
the Luttinger liquid model, associated to parabolic cor- 
rections to the linearized spectrum around the Fermi level 
17 24 . These effects lead to damping of the collective 



modes, which becomes relevant after a typical timescale 
idamp- For s-wave scattering in the weak coupling limit 
rj — /2'kvf ^ 1, idamp has been estimated to be 



id. 



amp 



512i?p/7r(r277)'^, with the Fermi energy. 



24 
Our 

approximation of neglecting such corrections is thus valid 




Fig. 3: Momentum distribution for right movers nk.ii(t) in the 
long-time regime (see eq. 171 as a function of momentum k. 
The parameters are the same as in fig. [2^. The inset presents 
the momentum distribution at a finer wavenumber scale close 
to the Fermi level. 



for pumping times t < idamp, while the spin-charge sepa- 
ration is best detectable in the "short-time limit" t < tg. 
We have toAdamp = Tr{nr]/EF)'^\n{2/jr])/256j. For the 
parameters of fig. [2^, io/^damp — 0.004. The regime of 
optimal visibility of the shoulders in the momentum dis- 
tribution is thus well described by our "non-interacting" 
collective modes approximation. Indeed, this picture is 
suitable to describe the "long-time regime" Iq < t < tdamp 
of fig. [3] as well, before damping yields saturation of the 
modes occupation at i > tdamp 
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A final issue to be addressed in view of the experimen- 
tal realization of our proposal is the role of finite tem- 
peratures, where our analysis above applies as well. The 
only differences are: (i) The presence of a non-vanishing 
Bose distribution in eq. (13 1. This yields a thermal seed 



for the amplification on top of the pure quantum fluctu- 
ations at T = 0, and decreases the time needed for the 
formation of well-resolved steps in nk^rif)- (h) The un- 
perturbed momentum distribution n\ ^ in eq. ( 14 ) has 
to be replaced with the finite temperature one, involv- 
ing a thermal smearing of order T around the Fermi level 
(for a ^ 1) on top of that purely induced by interac- 
tions. As thermal smearing involves wavenumbers up to 
order T/wp around /cp, the shoulders in the final momen- 
tum distribution at short times are thus better visible if 
v^qa ~ il/2 > T (a = c, s), which can be guaranteed 
by choosing a sufficiently large at a given temperature. 
In order for the Luttinger treatment to be reliable, the 
amplified qa should however be smaller than , which re- 
stricts the best choice of Q, to the window 2T < < 4i?p . 
The current experimental efforts to reach regimes of very 
low-temperatures T <^ Tp with cold fermionic gases would 
then further improve the frequency range for the best vis- 
ibility of the spin-charge separation. 

Experimental realization. — Our proposal could be 
experimentally realized with an equal mixture of quantum 
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degenerate fermionic K atoms confined into ID cigars in 
tlie two hyperfine states \F,mp) — |9/2,— 9/2) = | \.) and 
|F, ttif) = |9/2, -7/2) = | t)- Here, F is the total angular 
momentum and mp its projection along the quantization 
axis. It is important to realize that the atomic density 
corresponds to the charge channel, while the two hyperfine 
states above correspond to the (pseudo-)spin channel. As- 
suming the cigar of length L = 0.1 mm to be homogeneous 
and containing N — 10^ atoms, we have fcp = 3 x 10^ m^^, 
which corresponds to a Fermi energy and temperature of 
order Ep/h = 7 kHz and Tp ~ 60 nK, respectively. For 
neutral atoms, we assume contact s-wave scattering for 
which Vq = vJLf = as required by Pauli principle, and 
V(^/vf = nkpauj±/ Ep [25]. Here, a is the 3D s-wave 
scattering length and uj± the frequency of the transverse 
confining lasers creating the cigars. Notice that Ep ^ uj± 
justifies the effective ID treatment of the cigars. Assum- 
ing a scattering length of the order of a = 10 nm and 
uj± = 40 kHz, we obtain V^/vp = 0.5, which corresponds 
to the parameters in fig. [2^. 

As our proposal is optimized in the regime T < ri/2 < 
2Ep, and assuming T/Tp ~ 0.2, this corresponds to pump- 
ing frequencies in the range 3 kHz ^ ^2 ^ 28 kHz. Mea- 
suring the fermionic momentum distribution by a time- 
of-flight experiment [2], one should thus obtain a clear 
signature of spin-charge separation for pumping times of 
the order of 100 ms (for 57 — 10 kHz and 7 = 0.1), as ex- 
emplified in fig. [2^. Notice that pumping for larger times 
would lead to a situation similar to the one depicted in 
fig. [3] where spin-charge separation is much less clearcut 
and where temperature effects are likely to smear out most 
signatures of shoulders. 

Conclusion. — In this work we have shown the pos- 
sibility of parametrically amplifying collective modes in a 
modulated ID fermionic many-body system. The amplifi- 
cation is crucially affected by fermionic interactions which 
are here exactly treated within the Luttinger liquid pic- 
ture. This opens the perspective of similar observations 
in systems of higher dimensionality as well. 

Our analysis shows that the amplification of charge 
and spin density waves of the Luttinger liquid results in 
clear steps in the fermionic momentum distribution. The 
wavenumber extension of the steps directly reveals the dif- 
ferent momenta of the excited charge and spin modes and 
thus offers a tool for the detection of spin-charge separa- 
tion. In parallel, we show that the best resolution of the 
steps is achieved by modulations of relatively short times 
and that they survive thermal effects for large enough 
modulation frequencies. 

Our proposal of detection of spin-charge separation is 
particularly suitable for systems of cold fermionic atoms 
in ID optical lattices with modulated intensity. The 
fermionic momentum distribution is indeed the standard 
quantity measured in time-of-flight experiments. We 
stress that for our proposal no additional experimental 
setup is required on top of the already present tunable 



lasers creating the optical lattice. 

* ^ ^ 

We benefited from fruitful discussions with T. Bran- 

DES, T. ESSLINGER, T. GlAMARCHI, L. GOREN, and 

A. Stern. CDG enjoyed the hospitahty of the Indian 
Institute of Technology Roorkee while this work was 
completed. 

Additional remark: During the completion of this work, 
we became aware of ref. [17] where similar effects have 
been investigated. 
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